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Using the Teukolsky and Sasaki-Nakamura formalisms for the perterbations around a Kerr black 

' ■ hole, we calculate the energy flux of gravitational waves induced by a spinning particle of mass and 

^ ' spin S moving in circular orbits near the equatorial plain of a rotating black hole of mass M(3> /i) 

^ , and spin Ma. The calculations are performed by using the recently developed post-Newtonian 

QQ ' expansion technique of the Teukolsky equation. To evaluate the source terms of perturbations 

' caused by a spinning particle, we used the equations of motion of a spinning particle derived by 

. Papapetrou and the energy momentum tensor of a spinning particle derived by Dixon. We present 

CNj ' the post-Newtonian formula of the gravitational wave luminosity up to the order {v/cf" beyond the 

' quadrupole formula including the linear order of particle spin. The results obtained in this paper 
will be an important guideline to the post-Newtonian calculation of the inspiral of two spinning 
compact objects. 

o : 

^ '. 



I. INTRODUCTION 



One of the most promising sources of gravitational waves for kilo-meter size laser-interferometric detectors such as 
.5^ [ LIGO VIRGO 1^ and future laser-interferometric detectors in space such as LISA |^ is the coalescing compact 
. binary of neutron stars and/or black holes. Since it is a highly general relativistic event, detection of gravitational 
^ waves from those binaries will bring us very fruitful information about relativistic astrophysical objects if we know 
Ci ' physics of the final phase of the coalescence. Thus, there have been continuous efforts made by many authors to 
understand this phase 

Recently, it has been recognized that detection of the signal from a binary in the inspiraling phase is particularly 
very important because it can tell us a variety of parameters of the binary, i.e., mass, spin, etc. Furthermore it 

may provide some knowledge about the cosmological parameters j^,^ . In order to extract these important parameters 
from the data, we need an accurate theoretical template of the wave form. Especially, the accumulated phase of the 
emitted gravitational waves is very sensitive to the binary parameters. Thus, the rate of change in the frequency of 
the orbital rotation due to the radiation reaction of gravitational waves must be evaluated accurately. 

The post-Newtonian expansion is the standard method to calculate the wave form of gravitational waves. Recently, 
the energy loss rate to 2PN order beyond the quadrupole formula has been derived by Blanchet et al. Q and to 2.5PN 
order by Blanchet [ p^ for the binary composed of non-spinning compact bodies. The leading order effect of spin, 
which appears at 1.5PN and 2PN orders, have been evaluated by Kidder, Will and Wiseman ||ll| and Kidder p^ ]. 

However, dificulties and complications increase exponentially as one goes to higher orders wrEn the standardpost- 
Newtonian calculation technique. Hence it will be very useful if we have different approaches to the higher order 
approximation and if we are able to provide some guiding principle prior to the standard post-Newtonian calculations. 
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The black hole perturbation approach is the only one known method that is independent of the standard post- 
Newtonian approach and that can handle (not all but some important portion of) higher order post-Newtonian effects 
in a relatively straightforward way. In this approach, we consider a particle orbiting a black hole and assume the mass 
of the black hole M is much larger than that of the particle ^. 



The black hole perturbation approach is based on the perturbation equation derived by Teukolsky |13|, which 
applies to a general rotating (Kerr) black hole. One of great advantages of this approach is that it takes full account 
of relativistic effects by construction and numerical methods can be easily implimented to treat very general orbits. 
Moreover, it has been shown that one can formulate an analytical post-Newtonian expansion scheme in this approach 
as well. Poisson first noticed this fact and calculated the energy flux in the case of a particle in circular orbits around 
a non- rotating black hole to 1.5PN order flSj. A technically important point was to deal with the Regge- Wheeler 
equation, which is equivalent to the TeukoIsKy equation but has a much nicer property than it. Along this line, 
Sasaki developed a systematic method to proceed to the higher orders in the case of a non-rotating black hole [pT| , 
and Tagoshi and Sasaki jl^ gave the analytical expression for the energy flux up to 4PN order. Meanwhile Poisson 
calculated the case of a rotating black hole to 1.5PN order but directly dealing with the Teukolsky equation. 
Hence it seemed formidable to go beyond this order. Then extending the method of ref. a better method was 
developed by Shibata, Sasaki, Tagoshi and Tanaka to treat the case of a rotating black hole and the energy 
flux up to 2.5PN order was calculated. This was made possible by using the Sasaki-Nakamura equation [H, which 
is a generalization of the Regge- Wheeler equation for a non-rotating black hole. Recently the calculationhas been 
extended to 4PN order by Tagoshi et al. |20|. 

In all of these previous papers, the smallmass particle was assumed to be spinless. However, apart from the interest 
in its own light, for the purpose of providing a better theoretical template or at least a better guideline for higher order 
post-Newtonian calculations, it is desirable to take into account the spin of the particle. To incorporate this effect, 
we must know the energy momentum tensor of a spinning particle as well as the equations of motion. Fortunately, 
we can find them in the literature. The equations of motion of a spinning par ticle were first derived by Papapetrou 



n[, and they were put into more elaborate form by Dixon [2^ and Wald |23|. In particular, Wald ||2^ clarified all 
the conserved quantities along the general particle trajectory. On the other hand, Dixon ||22| succeeded to give the 
general form of the energy momentum tensor of a particle with multipole moments, which ofcourse includes that of a 
spinning particle as a limit. Hence, by using them, we can calculate the wave form and the energy flux of gravitational 
waves by a spinning particle orbiting a rotating black hole. 

Here, a word of caution is appropriate. Usually, we regard the small mass particle to be a model of a neutron star 
or a black hole. However, if we regard the spinning small mass particle as a Kerr black hole with mass /j, and spin 
angular momentum fiS, it should have the definite multipole moments ((x nS^). Since we neglect the contribution of 
these higher multipole moments in this paper, the particle in our treatment is not an adequate model for a (rapidly 
rotating) Kerr black hole. To incorporate the contribution of all higher multipole moments to represent a Kerr black 
hole is a future issue. Here, we concentrate on the leading order effect due to the spin of the small mass particle. 

This paper is organized as follows. In section 2 we briefly review the Teukolsky formalism. In section 3 we discuss 
the equations of motion and the energy momentum tensor of a spinning particle. In section 4 we solve the equations of 
motion to the linear order of the amplitude of spin. There we obtain a family of 'circular' orbits which have vanishing 
radial velocity and stay close to the equatorial plane. In section 5 we evaluate gravitational waves from the spinning 
particle in circular orbits and give the formula for the energy loss rate to 2.5PN order. In section 6 we consider the 
problem of the radiation reaction. There we show that the assumption that the orbit remains circular under radiation 
reaction is consistent with the energy and angular momentum loss rates in the linear order of spin. Then we evaluate 
the rate of change in the orbital frequency under this assumption. In section 7, brief summary and discussion are in 
order. 

We use the units G = c — 1 and the metric signature (—,-!-,+,+). The round (square) brackets on the indices 
denote the (anti-) symmetrization, e.g.. 



II. TEUKOLSKY FORMALISM 

In this section we briefly review the Teukolsky formalism. For details, see e.g., Ref. and references cited therein. 
In the Teukolsky formalism, the wave form and the energy flux of gravitational waves are calculated from the fourth 
Newman- Penrose quantity p5[ , which is expanded as 




Here, -2S^^{9) is the spin weighted spheroidal harmonics normalized by 

r \^2SZm' sin Ode = 1, (2.2) 

JO 
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and its eigen value is A. Then Remuj obeys the Teukolsky equation 



d / 1 dRi 



-V{r)Re„,^^Te„,Ur), (2.3) 



dr\A dr 
and 

+ Ai(r - M)K , , , 

V{r) = — ^- '— + ^iujr + X, (2.4) 

where A = — 2Mr + a? and K = (r^ + a^)ijj — ma. The source term Timuj{r) is constructed from the energy 
momentum tensor of the particle, and its explicit form is given later. 

The solution of the Teukolsky equation at infinity (r — > oo) is expressed as 

^'-^^ - 1^ = ^ ' (2-5) 

where r+ — M + V Af ^ — q,2 denotes the radius of the event horizon and R^l^^ is the homogeneous solution which 
satisfies the ingoing-wave boundary condition at horizon, 

r>in ^ f -OfmojA e , r > OO, ^.'i 

where k = uj — ma/2A/r+ and r* is the tortoise coordinate defined by 

dr* 



dr A 

For definiteness, we fix the integration constant such that r* is given explicitly by 

dr* , 
——dr 
dr 



(2.7) 



2M r+ , r — r+ 2Af r_ , r — r_ 

+ In — In , (2.. 



r+ - r_ 2M r+ - r_ 2M 



where r± — M ± V Af ^ — . 

Thus in order to calculate gravitational waves emitted to infinity from a particle in circular orbits, we need to 
know the explicit form of the source term Ti„icj{r), which has support only at r = tq where is the orbital radius in 
the Boyer-Lindquist coordinate, the ingoing- wave Teukolsky function Rg^uji''') at r = rp, and its incident amplitude 
^IniLu infinity. We consider the expansion of these quantities in terms of a small parameter, v'^ = M/rg. Note 
that V is approximately equal to the orbital velocity, but not strictly equal to it in the case of a ^ or S ^ 0. A 
systematic expansion method to calculate these necessary quantities has been developed in refs. p9| , pQ |, by considering 
the Sasaki-Nakamura equation first and then transforming the result to Teukolsky equation. 

In addition to these, we need to expand the spheroidal harmonics and their eigenvalues in powers of acu. Since 
oj = 0{il) where O is the orbital angular velocity of the particle, we have auj = 0{Mu!) = 0{v^). Thus the expansion 
in powers of ato is also a part of the post-Newtonian expansion. Note also that the spin parameter of the black hole 
a does not have to be small but can be of order M. 

The expressions of these quantities required to calculate the energy loss rate up to 2.5PN order are already obtained 
in ref. ||l9|. We digest the results omitting all the derivations. The homogeneous solutions of the Teukolsky equation 
with the ingoing boundary condition for £ = 2,3,4 are 

^-^2muj — — + — 2 — ttt: z 



30 45 1260 420 45360 11340 

/ — ,4 t 1*4 IILq^ 4:1^ All I r O-L t e 

-\-e I —^ —mqz — — z H — ^^^^ -f ^^^^^ mqz — ^^^^ z 



^ i 3 i ^ mqz^ 41 ^ 277 i g 31 z g 7 mqz' 

VT5" ~ ( 

9 

9/2: « 

e 



15 60 " 60 45 3780 22680 " 3780 1620 

2 , ~2 ^2 2 ^2 ^2 f. j-v-i /7 ^3 



n / z I , q z mqz t ^ mqz'' « 9 ■? i 9 9 c.\ ,„ 

^1 \ ma z + - - z - I a z m a z 1 (2 9) 

V30 40 60 240 60 30 90 120 /' ^ ' ' 

(2.10) 
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where e := 2Mll}, z := Lur and q := a/M. The incident ampHtudes are 



BS. - - 4 + Ki - 7 - log2) + + 0(,»)}, (2^4) 



- " {l + 0(e)}. (2.16) 





4mu. - g^2 



III. SPINNING PARTICLE 



To give the source term of the Teukolsky equation, we need to solve the motion of the spinning particle and also to 
;ivc the expression of the energy momentum tensor. In this section we give the necessary expressions, following refs. 
12],||J26|. 

lecting the effect of the higher multipole moments, the equations of motion of a spinning particle are given by 



D ,., ^ 1 

R-SP-[r)^2p^^\TX\T), (3.1) 

where w^(r) = dz^(T)/dT, r is a parameter which is not necessarily the proper time of the particle, and, as we will 
see later, the vector p'^ (r) and the antisymmetric tensor S^^{t) represent the linear and spin angular momenta of the 
particle, respectively. Here D/dr denotes the covariant derivative along the particle trajectory. 

We do not have the evolution equation for v^{t) yet. In order to determine v^{t), we need to impose a supplementary 
condition which determines the center of mass of the particle [ p2[ , 

5^''(r)p,(T) -0. (3.2) 

Then one can show that Pf^p^ = const, and S^j^nS^^ — const, along the particle trajectory [p3| . Therefore we may set 

pf" = fiu'' , u^u'^^-l, 
S""" = e^%„pPS'' , = 0, 

where fi is the mass of the particle, is the specific linear momentum, and is the specific spin vector w ith S its 
magnitude. Note that if we use instead of 5^" in the equations of motion, the center of mass condition (|]^) will 
be replaced by the condition p^S"^ = (see section 4). 

Since the above equations of motion are invariant under reparametrization of the orbital parameter r, we can fix r 
to satisfy 

u''{T)v^iT) = -1. (3.4) 



Then, from Eqs. (3.1, (3.2) and (3.4), ^''(t) is given as 



v>^iT) - u^{t) ^\[^^^ + \R^^U^{T))S^Hr)S^''{r)^ SP^t)R,,„^{t)u''{t)S'^-{t). (3.5) 

With this equation, the equations of motion (|3.1|) completely determine the evolution of the orbit and the spin. 
As for the energy momentum tensor, Dixon |E3| gives it in terms of the Dirac delta-function on the tangent space at 
= z'^{t). For later convenience, in this paper we use an equivalent but alternative form of the energy momentum 
tensor, given in terms of the Dirac delta-function on the coordinate space p6t|: 



(5^("(., .).«(., .)^^^^fc^)|, (3.6) 
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where v"{x, r), p°'{x, r) and S"^{x, r) are bi-tensors which are spacetinie extensions of u'^(r), P^{t) and S^'^{t) which 
are defined only along the world line, a;^ — (r).[] To define u"(x,z(r)), p'^{x,z{t)) and S°'I^{x,z{t)) we introduce 
a bi-tensor g°'^{x, z) which satisfies 



Yi^r^{x,z{T))^5"^, 
lim V/3r^(x,z(r)) =0. 



(3.7) 



For the present purpose, further specification of is not necessary. Using this bi-tensor g"^(a:,z), we define 

p"(a;,r), w"(a;,r) and S'"''(x,t) as 



P^{x,T)=r^X^,z(T))p^{T), 
V^{x,T)=r,X^,z{T))v^{T), 
S^\x,T)^r,X^,z{T))f^{x,z(T))S^^{T). 



(3.8) 



It is easy to see that the divergence free condition of this energy momentum tensor gives the equations of motion 
(3J). Noting the relations, 



V^5"(2^,^(r))<5(^)(x,z(r))==0, 



' i'^^){x,z{T)) 



d_ / ^W (x,z(r)) ' 



(3.9) 



the divergence of Eq. (3^) becomes 



V0r"/'(x)= / drr^(x,z(r)) 



5(4) (^-^(r)) 



^p''(r) + ii?^_(z(r))z;'^(r)5--(r) 



dr^ A-g\{x,z{T))f,{x,z(T)) 



(5W(x-z(t))~ 



-5 



S''^''(t) -2p[^(T)v''l(T) 



(3.10) 



Since the first and second terms in the right-hand side must vanish separately, we obtain the equations of motion ( [3.]] ). 
In order to clarify the meaning of p^ and S^^ ^ we consider the volume integral of this energy momentum tensor such 

/ Si^T^^dHf/^, where we take the surface S(to) to be perpendicular to u°'{tq). It is convenient to introduce a 

scalar function t{x) which determines the surface S(to) by the equation t{x) — tq, and dr/dx^ — —up at a; = 2(to). 
Then we have 



j^^^^9^^T-^dY.p = J d^x^g^5{T{x)-r,)-gJ^T-P{x) 



= J dr'^ir'-To) 
-P^(ro). 



p'^ +p 



2 dr 



(3.11) 



where we used the center of mass condition and the equation of motion for S^'^ . We clearly see p^^ indeed represents 
the linear momentum of the particle. 

In order to clarify the meaning of S^'^ , following Dixon 122], we introduce the relative position vector 



.= -gt'''d^a{x,z), 



(3.12) 



where (t(x, z) is the geodetic interval between z and x defined by using the parametric form of a geodesic y{u) joining 
z — y{0) and x = y(l) as 



^In the rest of this section, we use fi,u,a, - ■ ■ as the tensor indices associated with the world line z(r) and a, /3, 7, • • • as those 
with a field point x, and suppress the coordinate indices of z(t) and x for notational simplicity. 
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Then noting the relations 



lim = 0, Km = (3.14) 

it is easy to see that 



(3.15) 



Now that the meaning of S^"^ is manifest. ^From the above equation, it is also easy to see that the center of mass 
condition (^^) is the generalization of the Newtonian counter part, 

(fxp{x)x' = 0, (3.16) 

where p is the matter density. 

Before closing this section, we mention several conserved quantities of the present system. We have already noted 
that PfiP'^ = — /i^ and 5^5'' = are constant along the particle trajectory on an arbitrary spacetime. There will be 
an additional conserved quantity if the spacetime admits a Killing vector field 

^(m;.) = 0. (3.17) 

Namely, the quantity 

:= P% - ^^^''Cp;., (3.18) 

is conserved along the particle trajectory [Q. It is easy to verify that is conserved by directly using the equations 
of motion. 



IV. CIRCULAR ORBITS 

Let us consider "circular" orbits in the Kerr spacetime with a fixed Boyer-Lindquist radial coordinate, r = rg. We 
consider a class of orbits that would stay on the equatorial plane if the particle were spinless. Hence we assume that 
9 :— 9 — 71/2 ^ 0{S/M) ^ 1. Under this assumption, we write down the equations of motion and solve them up 
to the linear order in S. In appendix we give a further analysis in the case in which the spin vector (see below) is 
parallel or anti-parallel to the rotation axis of the black hole. 

In order to find a solution representing a circular orbit, it is convenient to introduce the tetrad frame defined by 

e^= (O, 0, O), 

e\^(0, 0, Ve, o), 

e\ = (-^sin0, 0, 0, !:!-^sin0), (4.1) 

where S = + cos^ 9, and = (e°(, e"^, e°'g, e"^) for a = ^ 3. Hereafter, we use the Latin letters to denote the 
tetrad indices. 

For convenience, we introduce oji ^ ojq to represent the tetrad components of the spin coefficients near the equatorial 
plane: 
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UJ 



31 



■^02 



■^32 



33 



0^ =c^i + O(02), 



"30 



"00 



"30 



■^13 



_ - Mr 

- ^2^1/2 ' 



a 



"33 
1 _ 



^12 



^11 



^4 + 0(6*^), := ^, 



605 + 0(02), ^5 



aAi/2 



"32 



(4.2) 



where 



, c _ fi y c 



Since the foUowing relation holds for an arbitrary vector 

e\^r = -i-r -^bc'v'T, 

the tetrad components of Df^/dr along a circular orbit are given explicitly as 

D 



D 



f'dr 
D 



where A, B, C, D and E are defined by|^ 



D 



<5zF = f' i-Bf eof) + o{e% 



A 
B 
C 
D 

E 



= UJ4V^ + UJ5V^, 



(4.3) 



(4.4) 



and we have assumed that = and — 0{d). 

Now we rewrite the equations of motion changing the spin variable. We replace the spin tensor with the unit spin 
vector, C", which is defined by 



or equivalently by 



oab o o-fy c /-a 



,a „,bccd 
bed 



ah ^,c/-d 



(4.5) 



(4.6) 



where tabcd is the completely antisymmetric symbol with the convention of £0123 = 1. As noted in the previous section, 
if we use the spin vector as an independent variable, the center of mass condition is automatically satisfied, while it 
becomes necessary to impose another supplementary condition 



CUa - 0. 



(4.7) 



■'■The symbols A '-^ E used here to define the auxiliary variable are applicable only in this section, and not to be confused with 
quantities defined with the same symbols, such as E for energy, in the later sections. 



7 



Then the equations of motion reduce to 



a b c o 73a 

— — — V u — bn , 



where 



-j- = 6.,,%"^ - Su^C'Rb, (4.8) 

i?'' R*\,dv''uX' = T^R\cdv'S^'', (4.9) 

and = \Rabeff^^ cd IS the right dual of the Riemann tensor. It will be convenient to write explicitely the tetrad 

components of R*abcd- Since we only need i?*bc(j at 0(6'"), the non- vanishing components of R*abcd given by 

^ 2^f^l23 ^ ~-fi^0213 — ^0312 = -^^1203 — "-^1302 = ~2"^2301 " + 0{9^). (4-10) 

Although wc do not need them, we note that the following components arc not identically zero but are of 0{6). 

^12121 ^1313' ^1010' ^2323' ^2020' ^^"^ ^3030' 
A. Lowest Order in S 

We first solve the equations of motion for a circular orbit at r = ro at the lowest order in 5*. For notational 
simplicity, we omit the suffix of in the following. For the class of orbits we have assumed, we have — Q and 
= O(0).Then the non-trivial equations are 

—v^=Av^ + Bv^^Q, (4.11) 
ctr 



d ^2 ~ ^ 







dr dr 



c 



3 




-C^=0, _ CM - U B e . (4.12) 



The equation ( 4.11 ) determines the rotation velocity of the orbital motion. By setting x := v^/v'^, we obtain the 
equation 

uji + 2uj2X + ojsx'^ ^ 0, (4.13) 

which is solved to give 

±VMr-a 



(4.14) 

The upper (lower) sign corresponds to the case that is positive (negative). Then, with the aid of the normalization 
condition of four momentum, v^v^ = —1 + 0(5*^), we find 

Note that, in this case, the orbital angular frequency is given by a well known formula, 

±%/M 

0= , ^ ■ (4.16) 

r3/2 ± VMo 



On the other hand, the equations of spin ( 1.12| ) are solved to give 



a sin(0 + ci) + /3c2 

c' = -a, I cM =Cii I cos(0 + ci) I . (4.17) 

-/3sin(0 + ci) — ac2 




8 



where C^, ci and C2 are constants, and 



A 



- A2 



VB2 - ^2 



it;", 



(4.18) 



The supplementary condition u^Cq = requires that C2 = 0. The condition CaC° = 1 imphes C± + C'\\ — 1- Further 
since the origin of the time r can be chosen arbitrarily, we set ci = 0. Thus, we obtain 

Here, we should note that flp ^ ft in general if a 7^ or S* 7^ (see below). 







/ a sin \ 


(?) 




cos 






y — /3sin0 y 



(4.19) 



B. Next Order 

Having obtained the leading order solution with respect to S, we now turn to the equations of motion up to the 
linear order in S. We assume that the spin vector components are expressed in the same form as were in the leading 
order but consider corrections to the coefficients a, (3 and Hp of order S. As long as we are working only up to the 
linear order in S, Eq. ( |3.5| ) tells us that v"' can be identified with In order to write down the equations of motion 
up to the linear order in S, we need the explicit form of which can be evaluated by using the knowledge of the 
lowest order solution, as 

i?° = i?^ = o{e), 



M 



= S^w^^C^ + 0(6*), 



.M 



i?2 ^ 2,'-^y^v^c^ + o{e). 



(4.20) 



First we consider the orbital equations of motion. With the assumption that = and = 0{9), the non-trivial 
equations of the orbital motion are 



i,i = Av° + Bv^ - SR^ = 0, 
^2 = {Cv° +Dv^)~e- SR^. 



(4.21) 
(4.22) 



The first equation gives the rotation velocity as before, while the second equation determines the motion in the 

^-direction. 

Again using the variable x — v^/v'~', Eq. (4.21) is rewritten as 



LUl + 2uJ2X + LU3X +3- 



-X ^ 0, 



(4.23) 



where Si 



SC±- The solution of this equation is 



It 



3S±VM 
2r3/2 



0(3^). 



(4.24) 



Using the relations {^A^), it immediately gives w" and v^. From the definition of the tetrad, we have the following 
relations. 



dt 



sm f 



asm 0— 

dr dr 



dt , , 9 , dip 
dr dr 



(4.25) 



Thus, the orbital angular velocity observed at infinity is calculated to be 
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dt 



a + x^/A 



+ 0(02) =±- 



M 



•3/2 ± aVM 



1 - 



3S± ±VMr-a 
2 ± aV Mr 



+ oie^). 



In order to solve the second equation ( 4.22 ), we note that = ^/T,6 ~ r6 and 

M 1 + 2.x2 



0{S). 



Then we find that Eq. (|4.22| ) reduces to 



~ Ml + 2x2- S\\M 
rO = r- —0 — 3- 



1 



■ COS( 



where S\\ = SC\\ ■ This equation can be solved easily by setting 6 — Oq cos Recalling that il^ = M/r^ 
obtain 



70 = 



(4.26) 

(4.27) 

(4.28) 
0(5), we 

(4.29) 



Thus we see that the orbit will remain in the equatorial plane if = 0, but deviates from it if ^ 0. We note that 

M 1 + 2a;2 



there exists a degree of freedom to add a homogeneous solution of Eq. (4.28), whose frequency, U,g = 

V r'^ 1 — X' 

is different from Op and which corresponds to giving a small inclination angle to the orbit, indifferent to the spin. 

Here, for simplicity, we only consider the case when this homogeneous solution to 9 is zero. Schematically speaking, 
the orbits under consideration arc those with the total angular momentum J being parallel to the z-direction, which 
is sum of the orbital and spin angular momentum J ~ L + S (see Fig. 1). 

Next we consider the evolution of the spin vector. To the linear order in S, the equations to be solved are 

The third equation is written down explicitly as 

= — 0^11 K sin (/) cos (/), 

with 

K := aD — pC — flpUj^r. 

Thus we find that 



C2 = -a + ^cos20. 



(4.30) 

(4.31) 
(4.32) 

(4.33) 



Since the spin vector S*" is itself of 0{S) already, the effect of the second term is always unimportant as long as we 
neglect corrections of 0(5^) to the orbit. 

The remaining three equations determine a, /3 and f2p. Corrections to a and /3 of 0{S) are less interesting because 
they remain to be small however long the time passes. On the other hand, the correction to fip will cause a big effect 
after a sufficiently long lapse of time because it appears in the combination of flpT. The small phase correction will 
be accumulated to become large. Hence, we solve fip alone to the next leading order. Eliminating C° and from 
these three equations, we obtain 



Then after a straightforward calculation, we find 



AC -ED 



^2 li ±VM (2r2 - 3Mr + a^) + ar^/^ {M - r) 



r3/2 r2 - 3Mr ± 2ax/Mr 



(4.34) 



(4.35) 



As noted above, fip is different from il for S± ^ 0. The difference f2p — il gives the angular velocity of the precession 
of the spin vector (see Fig. 1). 
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V. GRAVITAIONAL WAVES AND ENERGY LOSS RATE 



We now proceed to the calculation of the source term in the Teukolsky equation and evaluate the gravitational wave 
flux. For this purpose, we must write down the expression of the energy momentum tensor of the spinning particle 
explicitly. We rewrite the tetrad components of the energy momentum tensor in the following way. 



dr 



- {s^i^v'^si^Hx - z(r))) I 

= .,^Jdr ^ _ I . (5.1) 

The last line is the definition of A"'' and 

The source term of the Teukolsky equation is expressed in terms of the components of the energy momentum tensor 
projected with respect to the complex null tetrad defined as 

ra^ = [r + iacos6)~^ d — [e^ + ^^-t) ■ 

(5.2) 

We adopt the notation such as T„„ := n^vyT^n to denote the tetrad components. Then the source term is given by 



where 



T,„^ =4 / df2d^p-5p-i(B^ + B^*)e-""'^+''"*^^, 

V 27r 



B' = 



(5.3) 



/pA2i_i[p-VV+(p-2p-2A-iT,,,„)], 



2\/2 



B'* = 



--p8pAV+[p-V+(p-2pT^^)] 



2V2 



(5.4) 



with 



(r — lacoat 
m 



sm f 
iK 



— auj sm ( 



j cot 6, 



(5.5) 



and Q denotes the complex conjugate of Q. 

As we will see shortly, the terms proportional to S\\ in the energy momentum tensor does not contribute to the 
energy and angular momentum fluxes at the linear order in 5*. In other words, the energy and angular momentum 
fluxes are the same for all the orbits having the same S±, irrespective of the value of Si\ . Thus, we ignore these terms 
in the following discussion. Further we recall that the particle can stay in the equatorial plane if 5|| — 0. Hence we 
fix 6* 7r/2 in the following calculations. 

Using the formula (p^) , we obtain the amplitude of gravitaional waves at infinity as 



y ynn I ymn i ymm 



(5.6) 
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where 



ynn 



ymn 



iv/27r 



ujB 



d {u — mil) 



Ann - i^B, 



d_ 

dr 



ujB 



imu. 



1 

)=7t/2 rA 

-1 



i\/T^ ^ , f dt 



dr J 
1 

/2 %/A 



d 2iK 4 



dr A 



R 



-S {uj — mfl) 



Qj.2 



r A I dr 



d_ 

dr 

2iK K"^ 

"a7 ~ A2 



R 



imcu 



(5.7) 



and 



and 



A — - 



^rnn 



B^- - 



1 1 



{1 — S± {{2uJi + uj3)x + UJ2)} ■ 



1 



A _ 



1 



4r ^ 1 — 
i 1 



X + a, — \fKx, 0, —;=x + 1 I , 
A ^/A ' 



472 1- 



{2a; — {uix^ — 4a;2a; — wa) } 



4 - - — — - 



4V2r 1 - x2 V %/A 
1 1 



-x^ + aa;, - yA(l + x^), 0, — + x ) , 



{x^ + 5± (tJ2(l + 2x^) + waa;) } , 



(5.8) 



sinfl 



auj sin + j cot ( 



The Lorentz factor dt/dr which appears in Eqs. (5.7) can be calculated from Eqs. (4.25) as 

dt I 



dr r^/l — x'^ 



(5.9) 



(5.10) 



When the orbit is quasi-periodic, the Fourier component of gravitational waves do not have a continuous spectrum 
but takes the form, 



(5.11) 



Then the time averaged energy flux is given by the formula |19|] , 



dE\ 
'dt/ 



\ ^ I Zirnwn I 



fdE 
2-^ \ IT 



(5.12) 



The the z-component of the angular momentum flux is also given by the similar formula, 



dJ, 
dt 



GW 



E 



dJz 
dt 



(5.13) 
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In the present case of circular orbits in the equatorial plane, the index n degenerates to the angular index m and aj„ 
is simply given by mfi (n = m). Hence we eliminate the index n in the following discussion. 

Here we mention the effect of non-zero . If we recall that all the terms which are proportional to S\\ have the time 
dependence of e^'^^v'^ ^ we see that they give the contribution to the side bands. That is to say, their contributions in 
Zimuj are all proportional to 5{uj — mil, ± flp). Then, since the energy and angular momentum fluxes are quadratic in 
ZemujnJ they are not affected by the presence of S'n as long as we are working only up to the linear order in S. 

In order to express the post-Newtonian corrections to the energy flux, we define rjemuj as 



(-) 




(dE\ 




-A 


\dt) 



N 



where {dE/dt)^ is the Newtonian quadrupole formula 



dE\ 



5r^ 



32 



(5.14) 



(5.15) 



We calculate rjem up to 2.5PN order, i.e., to 0{v^). The result is 



r]2±2 



1 



107 
'21^ 



19 



/4784 



1 2 
— V 

36 

1215 
"896" 
5 

63 



^±1 



m±3 



-v 



q 



-V 



1 



m±4: 



m±2 



8064 
1280 4 

5 4 

3969 ' 

a/M and s 



1 

^12 
1215 
TT2 

40 
^189 

1 

1512^ 



12 



/3645 



V1323 
17 J_ 

^504 ^ 16 
1215 10935 



2q' 
1 



9qs V 



q — -qs V + 



18 



428 

vr H 

21 

793 
' 9072' 



4216 2134 



189 " 
535 

" loos' 



63 



V 448 
20, 



a H s 

^ 63 



1 



112 



17 



896 



1 



4032 



9072" 8064 



(5.16) 



where q := a/M and s :— Sj_/AI. The rest of rjim are all of higher order. We should mention that if we regard the 
spinning particle as a model of a black hole or neutron star, S is of order /x. Therefore the corrections due to S are 
generally small compared with the S*- independent terms in the test particle limit n/M <C 1. 
Putting all together, we obtain to 2.5PN order. 



dE\ 
di I GW 



dE 
~dt 



N 



1247 2 

V 

336 

44711 
" 9072 



47r - 



73 25 

— q ; 

12^ 4 

71 



33 

16'^' + y'^^ 



8191 
672 ' 



3749 
336 ' 



2403. 
112 ' 



(5.17) 



Since v is defined in terms of the coordinate radius of the orbit, the expansion with respect to v does not have a 
clear gauge-invariant meaning. In particular, for the purpose of the comparison with the standard post-Newtonian 
calculations, it is better to write the result by means of the angular velocity observed at the infinity. Using the 
post-Newtonian expansion of Eq. (4.26), 



the above result can be rewritten as 

1247 



[l-(^-s + q^ v'+^-qsv' + 0{v')^^ 



(5.18) 



dE\ 
dt I GW 




1 



336 

44711 
' 9072 



{MVl) 



2/3 



47r - 



11 



{MVl) 



+ ^qs ] {Amf/^ 



8191^ 
"672"' 



59 13 A,,, ^../c. 



(5.19) 
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where 





\ 32 


(-) 




N 





^ {]Vmy°^^. (5.20) 

Since there is no side-band contribution in th e present case, the angular momentum flux is simply given by 
{dJz/dt)fjy^ — {dE/dt)fjy^. The result ( |3.19| ) is consistent with the one obtained by the standard post-Newtonian 
approach to the 2PN order in the limit /u/M 0. The s-dependent term of order (Mfi)^/'^ is the one which 

is newly obtamed here. 

VI. RADIATION REACTION 

In this section, we consider the effect of radiation reaction on the orbit by equating the gravitational energy and 
angular momentum fluxes with their loss rates of the system. 

A. Conserved Quantities 

Here we consider the conserved quantities which are the first integrals of the equations of motion. First we give 
two conserved quantities which follow from the Killing vectors of the Kerr spacetime. The timelike Killing vector is 
given by 

A n asin6' ^ ^ , 

^^ = y^^^ + ^'^^ (6.1) 



and its derivative is 



2Af (H - cos-^ 6*) , n AMar cos 6 . , ^^ 
= + ^2 V-^] ■ (6.2) 



The rotational Killing vector is given by 



Xy. = a sm^ d\\ - + ^ ^ e^^, (6.3) 



and its derivative is 



X.;. = ((r - M)S + 2Mr^) e\e\^ 

2a-\/A sin cos 6* n n 2rsin^?\/A , n 

^ e [M^ ^ e [M^ H 



Then following the discussion around Eq. (3.18), we can construct the conserved quantities describing the energy and 
the 2-component of the angular momentum from these Killing vectors. They are given by 

E 1 

[a r, a sin 61 o Af (r^ - cos^ 9) S'^° 2Mar cos 9 S^^ 

= </ -\ 1;-^ -I i L I 

77 ''''^^ " 2^^'^'^^^"' 



a siw 



C0S6 / 9 . 9 „ , / 9 9n2\ S*^^ 



(a2sin2 0A-(r2+a2)')— , (6.5) 
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Since terms such as cos and cosdS'^'^ are higher order in 0, we neglect them in the following discussion. 

For a spinless particle, there exists one more conserved quantity on the Kerr spacetime, known as the Carter 
constant |2^. It is associated with the Killing tensor K^^^, which satisfies K(^f^^.,j-^ = 0. However, for a spinning 
particle, no such conserved quantity has been known. Nevertheless, one can show that there exists an approximately 
conserved quantity which corresponds to the Carter constant for a spinless particle. It is constructed as follows. 

In addition to the Killing vectors, the Kerr spacetime has an antisymmetric Killing Yano tensor, 

- 2a cos 9 e\e\^ + 2r e\e\^ , (6.6) 

which satisfies 

/m^;.) = 0. (6.7) 
Note that from /^^, we can construct the Killing tensor K^i, as 

K^,, = /^.// = r2g^, + 2EZ(^n,) . (6.8) 

When there is a Killing Yano tensor, the system possesses a quantity whose r-derivative is of 0(5*^), hence conserved 
to the linear order in S. Introducing the totally antisymmetric tensor, 

6 { ^ eV e\e\^ + e\e^] , (6.9) 



the approximate conserved quantity is expressed by |2^] 

_asin0 0^13 ^ ^1^30 _ 2^3^10 ^ a cos 5^3} 

-i.[^{acose(u^S^° - u^^^o + 2u°S^^) - ru'^S'"] . (6.10) 
/i V s 

The quantity corresponding to the Carter constant is 

C := 2Q - {J^ - aEf . (6.11) 
For the case of circular orbits under consideration, we find 

C 



,2 



'2aSi_ . (6.12) 



B. Frequency Shift due to Radiation Reaction 



The orbit of a non-spinning particle is completely determined by the three constants of motion, E, Jz and C . As 
mentioned above the counterpart of C for a spinning particle also exists in the linear order of S. This means that 
we need to calculate the radiation reaction to these quantities to obtain the orbital evolution of a spinning particle. 
In both non-spinning and spinning cases, the radiation reaction to E and Jz can be evaluated by equating their loss 
rates with the corresponding fiuxes emitted by gravitational waves. However, we do not know how to determine the 
reaction to C from the asymptotic behavior of emitted gravitational waves at infinity. Furthermore, there are also 
spin degrees of freedom in the present case. It is not clear at all how to evaluate the backreaction to the spin. In 
order to fully understand the radiation reaction, it will be necessary to derive some regularized radiation reaction 
force which acts on the particle. However, it is beyond the scope of the present paper. 

Although the rigorous evaluation of the change rate of spin seems formidable, there exists an order-of-magnitude 
estimate by Apostolatos et al. |29) , in which they calculate a torque acted inside a spinning star due to the radiation 
reaction force. According to their estimate, the rate of change of the spin is 0{v^^) smaller than that of the orbital 
angular momentum. This means we can safely ignore the time variation of the spin if w ^ 1. Since we are interested 
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in the case t; ^ 1 in this paper, we may then assume that the circular orbit obtained in section 4 remains circular with 
smaller r but with the same values of Sj^ and 5|| under the radiation reaction. Note that this implies the approximate 
Carter constant, Eq. (6.12), is conserved under the radiation reaction. 

If this assumption is correct, we can obtain an evolution sequence of the orbits under the radiation reaction. Namely 
the orbit is quasi-circular and slowly spiralling in with constant S'_l and S\\. Although not sufficient, we can consider 
the necessary condition for this assumption to be true by examining the consistency with the energy and angular 
momentum loss rates evaluated in terms of their fluxes emitted by gravitational waves. Here the consistency means 
that the relation 



dE 



dt I 



GW 



5J, 



dJz 
dt 



(6.13) 



GW 



holds, where the variations 6E and SJz are taken keeping the orbit circular with fixed S_l and S\\. As noted before, 
since there are no side-band contributions to the ener gy an d angular momentum fluxes in the present case, we have 
{dE/dt)GW — ^{dJ z / dt) GW ■ Therefore the condition ( |6.13| ) is equivalent to 



6E 



(6.14) 



In the following, we prove that Eq. ( |6.14| ) indeed holds. Since the S\\ does not appear in the expressions for E, J^, 
{dE/dt)GW and {dJz/dt)GW to the linear order in S, we only have to examine the situation with — 0. 
We introduce the function. 



ViE, Jz,r, S^) := -g'^UtUt - 2g''^utu. 



1, 



(6.15) 



which is guaranteed to be non- negative and it becomes zero when the orbit has no radial or vertical motion, i.e., 
v"^ = = 0. If we set E and Jz to the values for a circular orbit with r = tq, we have V — Q aX r = and V > aX 
any other points near r — rg. Hence dV/dr\r=ro = 0. 

The momentum components —put and pLU^ are different from the conserved energy and angular momenta due to 
the presence of spin. From Eqs. ( |6.5| ), we see 



-^iUt = E- ^iS^s {u°iE, Jz).uHE, Jz),r) , 
IJU^ ^ Jz~ fiS±j (u°(£^, lz),u^{E, lz),r) , 



(6.16) 



where 



M 



(6.17) 



Inserting these expressions to Eq. ( |6.15| ), taking the variation of V while keeping the orbit circular with fixed spin, 
and recalling dV/dr\r=ra = 0, wc obtain the relation between SE and 5Jz as 



2u* 
0. 



6E - ^iSi 



'dE 



SE 



de 
dJz 



SJz 



- 2u^ 



SJz - f^Si 



dE 



SE 



dj 
dJz 



SJz 



This gives 



SE 

SJz 



1 



de 
dTz 



dj 
dJz 



By using the following relations, 



1 



1 



de 
dE^ 



a )E — aJz 



dE' 



(6.18) 



(6.19) 



{Jz - aE) , 



1 + a 
r 




A 

1 a 

—u 
A r 



—u 

r 



r 



(6.20) 
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which hold i n the lowest order in S, it is easy to verify that th e term s in the square parentheses in the right-hand side 
of Eq. ( |6.19D become of 0{S^) or higher at r = vq. Thus Eq. (6.14) is shown to hold to the linear order in S and our 
assumption of stability of the quasi-circular orbit is found to be consistent. 

Under the assumption that the orbit remains quasi-circular with fixed spin, we can evaluate the frequency shift due 
to the radiation reaction by 



In) 



dE\ 

'^^ I GW 



(6.21) 



The post-Newtonian expansion of this quantity is calculated to become 



n^dt 96 M ^ ' 



336^ ' 



+ 



/ 34103 
V 18144 



81 , 
16^ 



79 



qs {Mn) 



4/3 



4159 
672 ' 



31319 
1008 ' 



809^ 
"84"' 



(Mn) 



5/3 



(6.22) 



VII. SUMMARY AND DISCUSSION 



In this paper we have investigated the gravitational waves emitted by a spinning particle in circular orbits around a 
rotating black hole. First we have solved the equations of motion of a spinning particle in the Kerr spacetime, assuming 
the spin of the particle is small and the orbit is close to the equatorial plain. Applying the Teukolsky formalism of the 
black hole perturbation, we have then calculated the first order corrections due to spin to the energy flux up to 2.5PN 
order. The effect of spin is always small {'^ 0{fj,/M)) compared with the spin-independent contributions if we take the 
limit /x/M 0. However, the result will be a useful guideline to the standard post-Newtonian calculations because 
our approach is totally different from the standard post-Newtonian approach and gives the leading spin-dependent 
terms in the fi/M expansion. Up to 2PN order, our results are in complete agreement with the previous ones obtained 
by the standard post-Newtonian method. The spin-dependent term at 2.5PN order, which we have newly obtained, 
will be verified by the standard post-Newtonian approach in future. 

In this paper we have restricted our analysis to a class of orbits which are circular and stays in the equatorial plane 
if the spin vector is orthogonal to it. In this case, the assumption that the orbit remains in this class under radiation 
reaction has been found to be consistent with the energy and the angular momentum loss rates evaluated from the 
gravitational wave flux at infinity. Although we have not considered it here, it seems possible to incorporate the 
orbital inclination. However, if the present restriction is relaxed, the orbit may become too complicated and the same 
technique used here may not work well. As long as one considers a spinless particle, the orbit is parameterized by the 
three conserved quantities, i.e., the energy, the z-component of the angular momentum and the Carter constant. In 
that case, the particle will sweep through a certain region of the phase space restricted by these conserved quantities 
sufhciently fast compared with the time scale of the radiation reaction. Hence the quasi-periodicity will be a good 
approximation. On the other hand, for a spinning particle, there are not sufficient number of conserved quantities 
to confine the orbit to a restricted region of the phase space. This means that the orbit is not guaranteed to be 
quasi-periodic. When we try to make a better template to be used for interferometric gravitational wave detectors in 
future by taking account of the effect of spin, this point may cause a difficult problem. 
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APPENDIX A: 



In this appendix, we present an exact solution of the equations of motion in the case when the spin vector S*" is 
parallel or anti-parallel to the rotation axis of the Kerr black hole. We assume the orbit is circular and lies on the 
equatrial plane. Thus the only non-vanishing component of the spin vector is SC,'^ = —S± and we set = = 
exactly. Further we assume — ~ 0, which will be found to be consistent. 

With these assumptions, the only non- vanishing component of i?" is 

= -S^^ {2v°u' + v\°)) . (Al) 
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Then all but one of the equations of motion are trivially satisfied and the remaining non-trivial equation is 



M , 



(A2) 



Besides the equations of motion there are constraint equations to be satisfied. We list them below. 

1. ^From the time derivative of the center of mass condition S°'^Uh — 0, we have the relation between u° and w° as 



Qah 

u^-^v"- SRb. 



which is equivalent to Eq. (3^). This gives 



(A3) 

(A4) 
(A5) 



2. The mass conservation Waw" = — 1 gives 



(A6) 



3. The normalization condition Uav"" = — 1 of is 



(A7) 



We may consider S±_ and r as freely specifiable variables. Then t he v aria bles to be determined are u^, ^ vP and v? . 
However there are five equations to be satisfied, i.e., Eqs. (A2), (A4) ~ (A7). Hence one might think the system is 
overdetermined and inconsistent. Fortunately this is not the case because Eq. (|A7| ) is guaranteed to hold by Eqs. ( A4), 
(A5) and ( |A6| ), as can be seen by contracting Eq. (AS) with Ua- Thus our assumptions, in particular = = 0, 
turn out to be consistent. 

In order to solve the above set of equations, we introduce the new variables. 



„3 y3 

Xv — and Xu — 



In terms of Xv and a;„, Eq. (A2) is rewritten as 

— Mr + aVA{xv + Xu) + AxyXu + aAlVA{2xu + x^) — 0, 

where 



a := = — s . 

r r 



(A8) 



(A9) 



(AlO) 



Multiplying Eq. (A4) by v?' and Eq. (A5) by vP , equating the right-hand sides of them and dividing it by vPv'^ , we 
obtain 



— cr2(2x„ + Xy){u°f = x„ - —a'^{2xu + Xy){u^f'. 
r r 



Using Eq. (AE), one readily sees this reduces to 



M 



(All) 



(A12) 



Thus we have obtained the coupled equations (A£) and ( A12 ) for Xy and Xu- They are solved to give 

_ -{2ra + 3Mra + aMa^) ± V4A/r3 -|- UaMr^a + ISlVPr'^a'^ + QaNPra^ - SAPra* + ^da?M^(j^ 

Xy 

_ r-Ma^ 
^" " r + 2Ma2^"- 



2VA(r - Mcr2) 



(A13) 
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Then from Eqs. (A6) and (A7), u^, and are found as 



^0 ^ V^EH , ^ ^'-V^-^' . (A14) 

1 XyXii 1 XyXy^ 

We note that the terms inside the square root of the expression for Xy in Eq. (A13) is not positive definite. Thus 
circular orbits do not exist for very large values of a (= S±/r). However, they always exist for physically reasonable 
values of a, i.e., for a < M/r < 1. We also note that the + (— ) sign in front of the square root corresponds to a 
co-rotating (counter-rotating) orbit if we restrict the range of a to be non-negative, i.e., < a < M. On the other 
hand, if we extend the range of a to —M < a < M, the ± signs become redundant. Here we take the latter option 
and take the -I- sign. 

Then a matter of interest is the stability of these circular orbits. In Fig. 2, we show the contours of radii of the 
inner-most stable circular orbits on the (a, S'j_)-plane. The contour of r = 6M passes through the {a,S±) = (0,0), 
which is the well-known minimum radius for a spinless particle in the Schwarzchild spacetime. One readily notices that 
the minimum radius increases as a increases for a fixed S± and it is larger for larger S± . Another interesting feature is 
that the minimum radius approaches r = M in the limit a — > — Af irrespective of the values of S±. Although this latter 
feature can be explained only in a fully relativistic context, the main feature of the contours can be understood as a 
consequence of the spin-orbit coupling, which is the dominant effect in a mildly relativistic situation. It is repulsive 
when the spin and orbital angular momentum vectors are parallel and attractive when they are anti-parallel. Now, if 
the contribution of the particle's spin to the spin-orbit interaction could be neglected, the contours of the minimum 
radii would be parallel to the S'j.-axis, with increasing minimum radii for larger a. On the other hand, if the particle 
were another Kerr black hole with the same mass M, spins of the black hole and the particle would contribute to the 
spin-orbit interaction in an exactly symmetric way, and the contours would be straight lines at 45° downward in the 
right direction. In reality, neither the spin of the particle can be neglected nor its contribution is as large as that of 
the black hole. This approximately explains the feature of the contours. 
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FIGURE CAPTIONS 



Fig. 1 A schematic picture of the precession of orbit and spin vector, to the leading order in S. The vector J 
represents the total angular momentum of the particle. The vector L is orthogonal to the orbital plane and 
reduces to the orbital angular momentum in the Newtonian limit. In the relativistic case, however, these vectors 
should not be regarded as well-defined. 

Fig. 2 The contours of radii of the inner-most stable circular orbits on the (o, S'_L)-plane. The conserved angular 
momentum is assumed to be positive so that the orbits with a > are co-rotating with the black hole and 
those with a < are counter-rotating. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/gr-qc/9602038vl 
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